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We describe the possibility of using LISA's gravitational-wave observations to study, with high 
precision, the response of a massive central body (e.g. a black hole or a soliton star) to the tidal 
gravitational pull of an orbiting, compact, small-mass object (a white dwarf, neutron star, or small- 
mass black hole). Motivated by this LISA application, we use first-order perturbation theory to 
study tidal coupling for a special, idealized case: a Schwarzschild black hole of mass M, tidally 
perturbed by a "moon" with mass (i < M in a circular orbit at a radius 6 > Af with orbital 
angular velocity Q. We investigate the details of how the tidal deformation of the hole gives rise to 
an induced quadrupole moment J?ij in the hole's external gravitational field at large radii, including 
the vicinity of the moon. In the limit that the moon is static, we find, in Schwarzschild coordinates 
and Regge- Wheeler gauge, the surprising result that there is no induced quadrupole moment. We 
show that this conclusion is gauge dependent and that the static, induced quadrupole moment for a 
black hole is inherently ambiguous, and we contrast this with an earlier result of Suen, which gave, 
in a very different gauge, a nonzero static induced quadrupole moment with a sign opposite to what 
one would get for a fluid central body. For the orbiting moon and the central Schwarzschild hole, 
we find (in agreement with a recent result of Poisson) a time-varying induced quadrupole moment 
that is proportional to the time derivative of the moon 's tidal field, = (32/45) M 6 Sij and that 
therefore is out of phase with the tidal field by a spatial angle 7r/4 and by a temporal phase shift 
7r/2. This induced quadrupole moment produces a gravitational force on the moon that reduces 
its orbital energy and angular momentum at the same rate as the moon's tidal field sends energy 
and angular momentum into the hole's horizon. As a partial analog of a result derived long ago 
by Hartle for a spinning hole and a static distant companion, we show that the orbiting moon's 
tidal field induces a tidal bulge on the hole's horizon, and that the rate of change of the horizon 
shape (i.e. the horizon shear) leads the perturbing tidal field at the horizon by an angle 4Mfi. We 
prefer to avoid introducing an ingoing null geodesic, as Hartle did in his definition of the phase 
shift, because the moon is in the central body's near zone (b <C and thus should interact 

with the horizon instantaneously, not causally. We discuss the implications of these results for 
LISA's future observations of tidal coupling, including the inappropriateness of using the concepts 
of tidal polarizability and tidal lag or lead angle, for the massive central body, when discussing LISA's 
observations. 

PACS numbers: 04.70.-s, 97.60.Lf, 04.30.-w 



I. INTRODUCTION AND SUMMARY 

A. Motivations 

One of the primary scientific requirements for LISA 
(the Laser Interferometer Space Antenna) is to map, 
in exquisite detail, the spacctimc geometries of massive 
black holes (and, if they exist, other massive, compact 
bodies) by using the gravitational waves emitted by inspi- 
raling white dwarfs, neutron stars, and small-mass black 
holes. This emission process has come to be called "Ex- 
treme Mass Ratio Inspiral" (EMRI, pronounced emm- 
ree). The possibility of making such maps from EMRI 
waves was discussed by Thorne in the early 1990s (e.g., 
in Rcfs. 0,0). In 1995 Ryan laid the first detailed 
foundation for such mapping: he showed that, when the 
massive, central body is general relativistic, axisymmet- 
ric, and reflection-symmetric, and the orbiting object is 
in a near-equatorial, near-circular orbit in the vacuum re- 
gion surrounding the body, the full details of the central 
body's metric are encoded in (i) the phase evolution of 
the waves and also in (ii) the evolution of the frequencies 



(or phases) of wave modulation produced by orbital pre- 
cession. Phinney Q has given the name "bothrodesy" to 
the mapping of a black hole's metric via EMRI waves, 
and bothrodesy has been identified, by the LISA Inter- 
national Science Team (LIST), as one of the prime goals 
for LISA p|. The initial phase of scoping out LISA's 
data analysis challenges for EMRI waves is now under- 
way H0. 

Ryan's proof that the EMRI waves carry a map of the 
central body's metric ignored completely the influence of 
tidal coupling between the central body and the orbiting 
object. Finn and Thorne have shown that, for values 
of the body and object masses in the range relevant to 
LISA, the tidal coupling can have an influence as large 
as a few percent on the evolution of the waves' inspiral 
phase — a phase that should be measurable to a fraction 
of a cycle out of tens or hundreds of thousands of cycles. 
Thus, the influence of the tidal coupling may be mea- 
surable with rather high precision. Because, in Ryan's 
analysis, the map is encoded redundantly in the EMRI 
waves' inspiral phase and in their modulations, it is rea- 
sonable to hope that the tidal coupling will break that 
redundancy in such a manner as to permit extraction of 
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both the map and details of the central body's response 
to the tidal gravitational pull of the orbiting object 

Thome has argued that if we are to keep an open 
mind about the physical nature of the central body from 
the outset [e.g., if we are to allow for the possibility that 
it is a boson star (e.g. 0,^2) or a soliton star (e.g. 0) 
rather than a black hole] , then we must describe the tidal 
coupling in a manner that can encompass all possible 
types of central bodies — a body-independent manner. 

In the case of the earth and moon, the tidal coupling 
is normally described in terms of the rise and fall of the 
earth's surface or ocean's surface, and in terms of en- 
ergy dissipation in the earth's oceans. Noticeably differ- 
ent from this, the tidal coupling in the case of a black 
hole has always, until now, been described in terms of 
the influence of the orbiting object's gravitational field 
on the hole's horizon — the perturbation of the hori- 
zon's 2-mctric (e.g. [3 Il5|). or the conversion of the 
tidal field into gravitational radiation at the horizon by 
gravitational blue-shifting and the energy and ang ular 
momentum carried inward by those waves (e.g., [lfj). 

One tidal feature in common between a black hole, 
the earth, a boson or soliton star, and all other conceiv- 
able central bodies, is the body's tidally-induced multi- 
pole moments and multipolar gravitational fields. It is 
these induced fields, acting back on the orbiting object, 
that change the object's orbital energy and angular mo- 
mentum, secularly change its orbit, and thereby alter the 
emitted gravitational waves. For this reason, Thorne [Tcj 
has proposed that we adopt these induced multipole fields 
or moments as our body-independent description of tidal 
coupling when analyzing LISA data. 

As a first step in exploring Thome's proposal, we com- 
pute, in this paper, the tidally induced quadrupole mo- 
ment and its back reaction on the orbiting object, in the 
special case where the central body is a Schwarzschild 
black hole, and the object is a distant, circularly orbiting 
moon. 



Framework and results 



when the central body is a black hole. 

To explore this, we study a model problem where 
the moon is orbiting circularly around a massive 
Schwarzschild black hole of mass M (3> n) at large dis- 
tance b (3> M). We assume the separation b is large 
enough that there exists an intermediate region between 
the hole and moon where (i) gravity is weak so space is 
nearly flat; (ii) the moon's tidal field does not vary ap- 
preciably. This region is referred to as the black hole's 
local asymptotic rest frame (LARF) ^8|. Because the 
spacetime is nearly flat, one can write down the full tidal 
field in the LARF (in Cartesian coordinates) to linear 
order in each multipole moment as |l9| 
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Here J^A t and i?A £ are the ^'th internal and external mo- 
ments; they are symmetric and trace free (STF) in their 
tensor indices = a\...ai |2jj. The "internal moments" 
J*A e characterize the central body, while the "external 
moments" J2A e characterize the gravitational fields of 
distant sources that perturb the central body. In our 
problem, the tidal field Ef** is physically the same as 
the external quadrupole moment ; they differ only by 
a constant scaling factor, £V, = — (l/3)i?.y. The inter- 
nal quadrupole moment is induced by the applied tidal 
field and characterizes the tidal deformation of the cen- 
tral body. 

Equation Q is the gravitational analogy to the multi- 
pole expansion of an electromagnetic field. It will be suf- 
ficiently accurate for our purpose, since we shall compute 
the nonspherical parts of the Ricmann curvature tensor 
by solving the linearized Einstein field equations. It will 
be shown in Sec. [H] that only multipole moments with 
^ = 0,2 are relevant to our problem. Dropping all other 
terms in Eq. Q and contracting with the unit spatial 
vector yields 



2M 
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Consider a moon of mass /i orbiting around a massive 
central body at a large distance. When the central body 

is a planet [HI (see Sec. Ill of Ref. [l^ for a review), where we have identified J in Eq. (JTJ as the total mass of 
the external tidal field produced by the moon, denoted the black hole and substituted the external tidal field Eff 
by £ f 7 xt , raises a tide on the central body and induces for =2y . In the last term, represents the quadrupole 

moment induced on the black hole by the external tidal 

field. 

In Secs. lIll| and llV A*l of this paper we compute the in- 
duced quadrupole moment in Regge- Wheeler gauge, ob- 
taining 



a quadrupole moment J^'" d 



that is proportional to £°f l ■ 



The proportionality constant is the body's polarizability. 
Because of viscous dissipation, the induced quadrupole 
moment J 2f /" d will be slightly out of phase with Sf**; it 
will have a small phase lag with respect to the applied 
field. This phase lag is generally referred to as the tidal 
lag angle, and can be defined equivalently as the ratio 
of the tangential and radial component of the tidally- 
induced force acting back on the moon. One objective of 
this paper is to explore whether this type of characteri- 
zation via polarizability and lag angle is also reasonable 
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The same result was recently derived by Poisson from 
calculating the averaged rate of change of mass and an- 
gular momentum of the perturbed black hole [2lJ. Note 
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that ^'" d is proportional to the time derivative of 
(a time derivative caused by the moon's motion) and is 
therefore completely out of phase with the external tidal 
field (by 90 degrees in time and 45 degrees in space) . As 
we will show in Sec. [V] this out-of-phase induced mo- 
ment is gauge invariant and is responsible for the torque 
that changes the orbital energy and angular momentum. 
Thus it is also responsible for the tidally-induccd portion 
of the orbital evolution and the phase evolution of the 
gravitational waves. 

The piece of the induced quadrupole moment that is 
proportional to and in-phase with the applied tidal field is 
ambiguous (in a sense that we shall discuss in Sec. II V B|l : 
in Schwarzschild coordinates and Regge- Wheeler gauge, 
it vanishes. If there had been an unambiguous piece of 
J >r i 1 " d in phase with Sf^, then this in-phase piece would 
have defined a polarizability, and the ratio of out-of-phase 
piece to the in-phase piece would have been, in a certain 
well-defined sense, the small tidal lag angle. Thus, our 
result can be regarded as saying that both the polariz- 
ability and the lag angle of a black hole are ambiguous 
(in the sense discussed in Sec. HVBll . 

Although wc find that the tidal lag angle in the LARF, 
in the case of a Schwarzschild black hole, is ambiguous, 
we can still define and calculate an angular tidal shift 
on the horizon (as opposed to in the LARF or out at the 
moon). We study this horizon phase shift in Sec.0 Har- 
tle [lfj has calculated 1 the tidal lag angle for the problem 
of a bulge raised on slowly rotating hole's horizon by a 
stationary moon, and he has shown it to be negative: the 
horizon's tidal bulge leads the applied tidal field due to 
the horizon's teleological definition (i.e., a definition in 
terms of the future fate of null rays). As in Hartle's case, 
we can compare the phase of the shape of our nonrotat- 
ing horizon to our moving moon's position by mapping 
the moon to the horizon with an ingoing, zero-angular 
momentum, null geodesic. In Sec. we find that this 
prescription leads to a lead angle between the moon and 
the horizon 



map = -Affi + fife* 



(4) 



where f2 is the orbital angular frequency of the moon 
and b* is the moon's tortoise coordinate b* = b + 
2Mlog(6/2M- 1). 

For comparison, Hartle's result for the tidal lead 
angle in the case of the rotating hole and distant, sta- 
tionary moon in the equatorial plane, is (after correcting 
a sign error, as discussed in footnote EJ) 



null map 
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4^, (5) 



Here a is the hole's specific angular momentum, and Qh 
is the horizon angular velocity. The radius of the moon's 



1 We review the principal results of Hartle's investigation in 
Sec. EH] 



position b is sufficiently large that the moon is essentially 
stationary. Throughout this paper, we use the super- 
script "H" to indicate results corresponding to Hartle's 
system, i.e., to a system with a stationary moon and ro- 
tating horizon. Other results (without the subscript "H" ) 
correspond to our system of a distant moon, orbiting at 
frequency Q, which perturbs a Schwarzschild black hole). 

Our result differs from Hartle's JSJ — even though 
we initially expected that the tidal phase shift should de- 
pend only on the difference in angular velocities of the 
applied tidal field and the horizon generators, so the re- 
sults would be the same. The terms that differ arise from 
the particular choice to map the moon to the horizon us- 
ing an ingoing, zero-angular momentum null ray. 

We prefer an alternative definition of the tidal lead 
angle, one that is independent of &*; we prefer to define 
the tidal phase shift as the angle between the perturbing 
tidal field at the horizon and the shear (which is the rate 
of change of the shape) of the horizon |23 ■ This definition 
avoids introducing null connections between the moon 
(which, at radius b <C f2 _1 , is in the near zone) and the 
horizon. Using this definition, we find that the shear of 
the central hole leads the perturbing tidal field at the 
horizon by an angle 



^Horizon — ^Horizon — 4MO. 



(6) 



The tidal lead angle is the same whether one considers a 
stationary moon perturbing a rotating hole or an orbiting 
moon perturbing a non-rotating hole. 

The rest of paper is organized as follows. In Sec. [HI 
we decompose the applied tidal field in the LARF into 
a time-dependent part and a static part. In Sec. IIIII 
we analyze fully the time-dependent part and deduce 
the dynamical part of the induced quadrupole moment 
[Eq. H24|)] . In Sec. lIVI we solve for the static perturbation 
and discuss the ambiguity in defining the static part of 
the induced quadrupole moment. In Sec.0 we study the 
phase shift between the deformation of the horizon and 
the applied tidal and compare the result with the phase 
shift as defined by Hartle. A brief conclusion is made 
in Sec. IVII Throughout the paper, we use geometrized 
units with G = c = 1. 



II. PROBLEM SETUP 

We study small perturbations of a non-spinning black 
hole caused by an orbiting moon. The unperturbed back- 
ground metric is the Schwarzschild metric: 



ds 2 



1 



2M 



■) 



dr 2 



dt2 + -, K77T +r 2 (de 2 +sin 2 8dtf ), 

1 — Zl\l I v 

(7) 

where M is the mass of the central hole. At large radii 
(i.e., in the LARF), we will study the perturbations 
in a notation that treats the Schwarzschild coordinates 
(r, 6, (f>) as though they were flat-space spherical coordi- 
nates. These coordinates are related to the Cartesian 
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coordinates (x,y,z) = (x 1 , x 2 , x 3 ) by 

(x , x 2 , x 3 ) = r(sin 9 cos <f), sin 9 smtft, cos 9). 

We will denote the radial vector with length r by x, the 
unit radial vector by n, and their components by x 3 and 
n J , respectively. 

Let a moon of mass [i move along a circular orbit with 
radius b in the equatorial plane (b ^S> M 3> fx). The 
moon's position is specified by 



x s = b n s = b (cosOt, sin fit, 0), 



(8) 



where the superscript "s " stan ds for the "source" of the 
perturbation and fl = y/M/b 3 is the moon's orbital an- 
gular frequency, satisfying Qb -C 1. The moon's tidal 
field is the double gradient of the moon's Newto- 
nian gravitational potential. Its value in the LARF (at 
r <C b but r>M) is well approximated by 



ccxt 
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1*J 
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r=0 



(9) 



Note that although the applied tidal field is defined in 
the LARF, the induced quadrupolar field J^lf of greatest 
interest is not in the LARF. but further out in the vicinity 
of the moon's orbit, where it interacts with the moon. 

The tidal field El can be decomposed into spherical, 
harmonic modes |23j. The result of the decomposition is 



the outer boundary condition for the problem: the O(r ) 
terms in the perturbed tidal field £^ must go to £?** in 
the LARF [Eq. ©]. 

The inner boundary condition is set differently, de- 
pending on whether the perturbations are static or time- 
dependent. For the static perturbations generated by 
£if, we impose a "regularity boundary condition": the 
perturbations must be physically finite at r = 2M . For 
the time-dependent perturbations generated by £^' 2 2 , 
we impose the "ingoing-wave boundary condition" : the 
perturbations have the asymptotic behavior ~ e =F4wr 
when approaching the horizon. Here r* is the tortoise 
coordinate r* = r + 2Mlog(r/2M - 1). 



III. TIME-DEPENDENT PART OF THE 
PERTURBATION 

A. The perturbed metric 

We will specialize to (£, to) = (2, 2) in solving for the 
time-dependent part of the metric perturbation. The 
(I, to) = (2, —2) results can be obtained by complex con- 
jugating the (2, 2) results. For briefness, a superscript 
"22" will not be added to quantities calculated in this 
harmonic mode in this section, unless a distinction is 
needed. Throughout this section, we refer to AppendixlBl 
for details of the perturbation calculation. 

In the standard Rcggc- Wheeler gauge, the (£, to) = 
(2, 2) time-dependent perturbations take the form |24| 



£ff = 




O1/-20 _ 0fs22 -iu>t 



ooxt.20 , ccxt,22 



(10) 



with uj = 2f2 and £^ ' (to = 0, ±2) equal to the 
corresponding < 3( i 2rn term. Here the ( 3^ 2 f a are position- 
independent, rank-2, symmetric trace-free (STF) tensors 
defined in Eqs. (|A3|I - (|A4|) and are related to the familiar 
t = 2 spherical harmonics Y~ 2m (6>,0) by Eq. (See 
Eqs. (2.7)-(2.14) of Ref. [2^ for the general mapping be- 
tween order I spherical harmonics and rank-£ STF ten- 
sors). The explicit values of the tidal field components 
are 
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r 2 K sin 2 9 
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Here H, Hi and K are functions of r alone. These radial 
functions are solutions of the perturbed Einstein equa- 
tions; they can be constructed from the Zerilli function 
Z(r) [2f|, which satisfies a second order ordinary differ- 
ential equation [Eq. i|B4|) ]. Specifically, Hi, K and H 
are given in terms of Z(r) by Eqs. I|B1|) - HB3|) . Instead 
of solving for Z(r) directly, one may obtain the Zerilli 
function from its odd-parity correspondent, the Regge- 
Wheeler function X(r), which obeys a simpler differential 
equation 0iEl| that ls easier to solve [Eq. 1B5|) ]: 



ar \ r J V r r a / 



X(r) = 0, 



The tidal field £ff l [Eq. fTOJ] is the source of pertur- 
bations of the central hole; it is an even-parity £ = 2 
external tidal field. We shall therefore perform our cal- 
culation in the even-parity Regge- Wheeler gauge, mode 
by mode (£ = 2, to. = 0, ±2). The tidal field £ff also sets 



where d/dr* = (l-2M/r)d/dr. The Zerilli function Z(r) 
is expressed in terms of X(r) by Eq. I|B6(I . Thus, the 
metric perturbation is determined by the single radial 
function X(r), by way of Eq. (|B6(1 to get Z(r) and then 
Eqs. l|BT )l -l|n5 ]l to get H x , K, and H. 



5 



The analytic solution for X(r) with the ingoing-wave 
boundary condition at horizon was derived by Poisson 
and Sasaki [2]}. Their solution. X H in their notation 
and for the limiting case wr< 1, is what we have used in 
our analysis. With our slow motion assumption Qb <C 1, 
X H (ojr <C 1) will be sufficient to cover the region inside 
the moon's orbit — including the LARF, where we read 
out the induced quadrupolc moment. Following Poisson 
and Sasaki's notation, we define the dimcnsionless quan- 
tity 



2Mw. 



(13) 



We then combine Eqs. (3.4), (3.11), and (3.12) of Ref. |2J 
to obtain 



X>r«l) = A(^) , 



x F(ci,c 2 ;c 3 ;l 



iu(r-2M) 

r 
2M 



(14) 



where A is an overall scaling factor that did not appear 
in Ref. |27j but will be determined by the outer boundary 
condition in our problem; F is the hypergeometric func- 
tion with parameters [Eq. (3.11) of Ref. (23 with I = 2] 



ci = — ie + 0(e 2 ), 
c> 2 = 5-ie + 0{e 2 ), 
c 3 = 1 — 2ie. 



(15) 



Note that expression (|14fl for X H is only accurate to first 
order in e. We then expand Eq. 1|14[> in large r and keep 
terms to first order in e 
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where f = r/2M is the dimcnsionless radius. Next, we 
use Eq. (|B6J) to get Z(r). Then the perturbed metric com- 
ponents can be constructed using Eqs. (|B1|) - (|B3(I . In the 
following all quantities will be calculated up to first order 
in £ and we will suppress "0(e 2 )" in our expressions. 



From this equation we can read off £. 



(1) d(1) fi p c rs + 

order tidal field in Cartesian coordinates in the LARF, 
from the Riemann tensor in Schwarzschild coordinates. 
By the procedure outlined in this paragraph we have de- 
duced the following (£,m) = (2,2) part of £y in the 
LARF: 



f(l).22 
■'ij 



SA 

' 4M 3 



-te - 



■E 

n— 5 



IE 



(19) 



The outer boundary condition states that the O(r ) 
[i.e. O(r )} term of S^' 22 must equal S'rf- 22 [Eq. (tTlbl l. 
This determines the scaling factor to be 



(20) 




Inserting Eq. (|20|l into Eq. (|19|) . we can write £. 



;i),22 
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c,(l),22 _ £,ext,22 
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(21) 



Here the 0(l/f 5 ) term, by Eq. (J2J, contains the in- 
duced quadrupole moment. The 0(l/f 6 ) and higher 
terms are proportional to the 0(l/f 5 ) term and con- 
tain no new information; they represent the non-linear 
coupling between the induced quadrupolc and the black 
hole's monopolc moment. 

Comparing the 0(1 /f 5 ) term in Eq. @ and the 
0(l/r 5 ) term in Eq. (J3TJ, we find that 



_ —M 6 f cxt ' 22 
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Complex conjugating this equation yields the (£, m) 
(2, —2) part of the induced quadrupole moment: 



2-2 
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(23) 



Thus, the time dependent part, i.e. the dynamical part 
(DP), of the induced quadrupole moment is given by 



B. Induced quadrupole moment in the LARF 

Once the perturbed metric is known, it is straightfor- 
ward to calculate the full Riemann tensor and extract 
from it the first-order tidal field in the LARF: 
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where a superscript of "(0)" or "(1)" indicates the quan- 
tity is of zeroth or first order in the perturbation. In 
our calculation, we found it convenient to look at the 
OrOr component of the first-order Riemann tensor in the 
LARF, since 



R {1) -R {1) n n- 
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(24) 

This ag rees with the result recently obtained by Pois- 
son |2J by a very different method. Note that the in- 
duced quadrupole moment is proportional to the time 
derivative of the applied tidal field. Hence the induced 
quadrupole moment and the applied tidal field are com- 
pletely out of phase with each other (ir/4 phase shift in 
space, it/ 2 in time). This leads to a dissipative force 
acting back on the moon. 

From the induced quadrupole moment Q24[l. we define 
a corresponding Newtonian potential in the LARF and 
out to the moon's orbit: 



S^ind, dp n l n 3 



(25) 
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Then the force acting back on the moon can be found by 
evaluating the gradient of <I> at the moon's position: 



F = -V*L. = - 



32 //x\2 / M 
"5\bJ \T 



13/2 



(26) 



Equation (|26() shows that the force is tangential and op- 
posite to the moon's motion. The energy loss from the 
moon's orbital motion is then 



E 



•32 



(27) 



where v = Qbe^. It is straightforward to show that 
there is also an angular momentum loss of magnitude 
E/fl. Equation l(2T|) agrees with Poisson and Sasaki's 
calculation of the rate at which the perturbation carries 
energy into the black hole's horizon at the leading post- 
Newtonian order 1271. 



IV. THE STATIC, AXISYMMETRIC PART OF 
THE PERTURBATION 

A. Static induced quadrupole moment 

We now specialize to the even-parity, static part of 
the moon's perturbation: (£,m) = (2,0). The Reggc- 
Wheeler metric for this type of perturbation has the 
form l2l 



,(20) 
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where "Diag" is short for diagonal matrix and H 2 and 
K 2 are functions of r only. The general solution to the 
field equation governing H 2 can be exp ressed in terms of 
the associated Legendre functions [lj| : 



H 2 (r)=a 2 Pr 
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(29) 



where a 2 and (3 2 are constants to be determined. The 
solution to K 2 can then be obtained from that of H 2 
(Appendix C). As r approaches 2M, we have 

Qe(r/M-1) ~ (r/A/r 1/2 , 

so the Q 2 term in Eq. (|29|l becomes divergent at r — 2M 
and we must set the coefficient (3 2 to be zero in order for 
the perturbation to be finite there. As r goes to infinity 2 , 



P?(r/M-l)~(r/MY 



Valid for all Re I > -1/2. 



Therefore the remaining P 2 term in Eq. I|29|) keeps grow- 
ing quadratically as r becomes large, corresponding to 
the non-asymptotic-flatness due to the presence of the 
moon. 

With the metric perturbation , we compute the 
Riemann tensor from the full metric and series expand 
the result up to linear order in a 2 (i.e., first order in 
the perturbation) . The OrOr component of the resulting 
first-order Riemann tensor is found to be 



(i) 

OrOr 



3«2 



Y 



20/ 



(30) 



From this and from Eq. i|10|) . we obtain the first-order 
tidal field in the Cartesian basis 



7(1) 



3a2 0^20 
M 2 11 ' 



(31) 



The static, first-order tidal field thus contains only an 
O(r ) term, which should be identified as the static part 
of the applied external field £ ' [Eq. (| 1 1 ap ] . The co- 
efficient a 2 is determined from this identification to be 
Q'2 = 1/471-/45 fiM 2 /b 3 . Since there is no 0(l/r 5 ) term 
present in Eq. (|31|l . we infer that there is no static in- 
duced quadrupole moment: 



(32) 



This is quite a counter-intuitive result. It is worth 
pointing out, however, that the absence of negative pow- 
ers of r in Eq. H31JI follows directly from the regular- 
ity condition we imposed at r = 2M. If the radius 
r = 2M were well inside the central body itself, which 
naturally is the case for any nonrclativistic body with 
weak self-gravity, then the Q 2 term in Eq. (|29(l would 
survive and give rise to an induced quadrupole moment. 
Equation (|32|l may also be the consequence of the gauge 
(Rcggc- Wheeler) we choose to work in. Is it possible to 
give a gauge-invariant definition of static, induced mul- 
tipolc moment in a non-asymptotically-flat spacetimc? 
This is the question we shall investigate in the next sub- 
section. 

We now summarize and conclude that the total in- 
duced quadrupole moment in our chosen gauge is 



ind 



20 



22 



1-1 



45 M & « ' 



(33) 



which is proportional to the time derivative of the exter- 
nal tidal field — not the field itself as one would expect 
for Newtonian tidal couplings. 

Lastly, we move from the LARF to the perturbed hori- 
zon and examine the effect of the static perturbation 
there. Hartle has shown that to first order in the 
perturbation, the coordinate location of the event hori- 
zon of a slowly rotating black hole perturbed by a sta- 
tionary distribution of matter is still at r = 2AI. This 
is also true for a Schwarzschild black hole under static 
perturbations. Evaluating the full metric at r — 2M, we 
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find the horizon metric is given by 



changes the hole's momentum 3 : 



ds 2 H 



4AP 



1 - 2^(M 2 /b 3 )P 2 (coa9) 



(d6 2 + sin 2 9d(j) 2 ) , 
(34) 

where P 2 is the Legendre function. From this metric the 
scalar curvature of the horizon is obtained as 



f 



2AP 



1 - 4/i(M 2 /& 3 )P 2 (cos6>) 



(35) 



So it is clear that the shape of the horizon does acquire 
a small quadrupolar component. But this deformation is 
not accompanied by an induced quadrupole moment in 
the LARF, at least in our chosen gauge. 



B. Ambiguity of the static induced quadrupole 
moment 

In the previous subsection, we found that a 
Schwarzschild black hole has a vanishing static induced 
quadrupole moment (SIQM) in response to the exter- 
nal tidal field £° xt ' 20 . To see that this vanishing of 
the SIQM might possibly be a gauge effect, imagine re- 
placing the radial coordinate r in the expression $ = 
(l/2)£^' rfn^r 2 for the external tidal Newtonian po- 
tential by r = f(l + xM 5 /^ 5 ) 1 / 2 , where \ lS some di- 
mcnsionless number of order unity. The result is $ = 
(l/2)£!f ' 2t W'f 2 + (x/2)M 5 £^ xt ' 2 Vn7f 3 . By com- 
paring this expression with Eq. Q we read off a SIQM 
= (x/18)M 5 f J ° J xt ' . In Newtonian theory this pro- 
cedure would obviously be naive, but in general relativ- 
ity, where the unperturbed black hole metric can be ex- 
panded in powers of M/r and the coefficients in that ex- 
pansion depend on one's choice of radial coordinate and 
that choice is a "gauge" issue, this type of procedure is 
not obviously naive at all. 

From our point of view, the best way to explore the 
gauge dependence of the SIQM is to ask whether it is 
physically measurable. If (as we shall find) physical mea- 
surements give a result that is ambiguous at some level, 
then that ambiguity constitutes a sort of gauge depen- 
dence of the SIQM. 

In this section we shall study a thought experiment for 
measuring the SIQM, one based on coupling to a small, 
static external "test" octupole field £ijk (proportional 
to the symmetrized and trace-removed gradient of some 
fiducial external quadrupolar tidal field). For simplicity 
we take £ijk to be axisymmetric around the same z axis 
as our static external tidal field £, cxt ' 20 ; i.e. we take it to 
be proportional to a tensor spherical harmonic of order 
(e,m) = (3,0): 



£ 



ijk 



9®- 



The analysis in Ref. [2!j says that, any SIQM J\j (cre- 



ated in the black hole by £ ? x ' ) will couple to the exter- 
nal octupole moment to produce a force that gradually 



(36) 



(Eq. (1.12) of Ref. [3; Eq. (4b) of Ref. j^). The same 
will be true if the central black hole is replaced by a 
neutron star or any other spherical body. The rate of 
change of momentum P l can also be evaluated by a sur- 
face integral of the Landau-Lifshitz pseudotensor i^L m 
the LARF fiji: 



P' 



(-3) *LL dS 3 



(37) 



Equations (|36|l and l|37|) for the coupling-induced force 
on the hole actually have ambiguities that arise from 
nonlincarities in the Einstein field equations. The ori- 
gin of those ambiguities is discussed with care in Sec. I 
of Thorne and Hartle 01- I n this subsection we use 
Eq. Ij37(l to calculate the force on the hole, and shall 
identify the ambiguities as those terms in which the force 
depends on the location of the integration surface. The 
result of our calculation will tell us, by comparison with 
Eq. 13tj|) . the SIQM and the amount by which it is am- 
biguous. 

To compute the pseudotensor for insertion into 
Eq. (|37f) . we must solve for the metric perturbation con- 
taining both the quadrupole and octupole terms: 



,(20) 



,(30) 



Diag[(l-2M/r)J3i, - 



Hp 



e=2,3 



2M/r ' 



(38) 



r 2 Ki, r 2 Kism 2 xF w ((9,0). 

When £ = 2, the general solution to H2 is given in 
Eq. For t = 3, we have 



^)= a3 Pi(^-l) + P 3 Qi(^-l) 



(39) 



In order that both types of perturbation be finite at r = 
2M, fa and fa must be set to zero [see Eq. I|30|l ]. In order 
to deal with more general cases, however, we keep non- 
vanishing values for fa and fa in Eqs. (j29(l and H39[) so the 
following analysis will be valid for stars as well as black 
holes. [For central bodies other than black holes, fa and 
«2 (and similarly fa and 0:3) are not independent of each 
other: fa is proportional to «2 with a proportionality 
constant that depends on the body's internal physical 
properties]. Having specified the metric perturbation, 



3 The test octupole field may also induce a static octupole moment 
J?ijk m the central black hole, which will couple to the external 
quadrupolar tidal field. This coupling, as we shall show, con- 
tributes to the ambiguities in the definition of the SIQM. 
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we then insert the full metric into the expression for the 
pseudotensor (Eq. (20.22) of Ref. [H) 



16tt l J 



,A0" 



( 



x (2g vp g aT -9 P ,gu T )s VT ,xS'" T ^}, (40) 

and evaluate the surface integral at some radius r = R in 
the LARF. Because of the axi-symmetry of the perturbed 
spacetime, only the ^-component of P 1 is nonzero. The 
result, up to first-order coupling and with uninteresting 
numerical coefficients being suppressed, has the following 
form: 



P z 



a 3 a 2 



cczfa 



M 4 M 3 



— & — 

M 2 M 



1 & 



(41) 



faa 2 



M \ M \ 



where "&" is to be read "and a term of the order" . 

The constant terms in Eq. I|41|l [i.e., the "l"s] that are 
independent of the integration radius R are the ones to be 
compared with Zhang's result l|36|) so as to deduced the 
gauge-invariant SIQM. Other terms that depend on R 
constitute ambiguities 4 . Terms with positive power(s) of 
R/M appear because the spacetime is not asymptotically 
flat, and they prevent us from minimizing the ambiguities 
by simply pushing the integration surface to infinity. 

Let us step back and write down the most general form 
that the SIQM can take. By order of magnitude analysis 
of the response of any physical body (black hole, star, 
planet, moon, ...) to a tidal field, we must have 



a. 



L 5 (l+0£ 



cxt,20 



(42) 



Here L is the size of the body (L ~ M for a black hole) 
and £ is a dimcnsionless number describing the SIQM's 
dependence on the integration radius R — deviations 
from being well-defined. From Eq. I).'U[I. we know the 
external tidal field scales as ~ a 2 /M 2 . Similarly for the 
external octupole field, £ijk ~ a 3 /M 3 . Using these rela- 
tions, Eq. (|36fl becomes 



L 5 (l 



a 2 a 3 



L 5 £L 5 ~ 
M5 + M5 



(43) 



4 The M 2 /R 2 term includes the effect of any octupole moment 
induced on the central body by the test octupole field. Note 
that P l is a dimcnsionless vector. On dimensional grounds, the 
coupling between any induced octupole moment and the external 
tidal field must take the form J rvk £ji e /R !l to contribute to P 1 . 
(Nonlinear coupling to the monopole moment can lead to similar 
terms that scale as higher, but not lower, powers of 1/R). 



Here the first term in the square bracket should be identi- 
fied as the "l"s in Eq. (|41|l [note again that fa an d ol 2 are 
not independent of each other for stars] ; and the second 
term should be identified as the sum of all independent 
terms: 



M 5 



d4 d A/T 

^-j & ... & & — & .... (44) 
M 4 MR v ' 



In the case of a black hole we have L ~ M and the 
smallest the right hand side of Eq. (|44|) can be is ~ 1 (for 
R ~ M), so £ > 1, i.e. the SIQM for a Schwarzs child 
black hole is ambiguous by an amount > M 5 £°* t ' , i.e. 
totally ambiguous, since the largest we could expect 
to be is~M 5 £f' 20 . 

For central objects with L 3> M (e.g., the Earth) we 
must choose R > L. The right hand side of Eq. I|44|) is 
then minimized by setting R ~ L, giving ^ ~ M/L <C 
1 (~ in the case of the Earth) for the fractional 

ambiguities in the SIQM. 

We comment that our result for a Schwarzschild black 
hole differs from what Suen has derived. Suen has given 
an unambiguous prescription to read out static multipole 
moments in non-asymptotically-flat spacetimes, which is 
based on transforming coordinates into a particular set 
of de Donder coordinates [3l||. He has used his prescrip- 
tion to calculate the induced quadrupole moment of a 
Schwarzschild black hole when it is perturbed by a static, 
equatorial matter ring at large distances 01- According 
to his prescription, the SIQM does not vanish. It is pro- 
portional to the tidal field produced by the ring: 



a. 



•— M 5 £ Tias 
21 ^ 



(45) 



The incompatibility between this result and the van- 
ishing SIQM that we derived in Sec. IIV Al in Reggc- 
Whcclcr gauge and Schwarzschild coordinates illustrates 
the ambiguities of the SIQM. Both results, zero and 
(— 4/21)A J f 5 £jj are less than or of order the ambiguity. 



V. THE TIDAL PHASE SHIFT 

In the LARF, the time-dependent induced quadrupole 
moment is 7r/4 out of phase with the perturbing tidal field 
f Sec. 1111 A"|) . This large phase shift is quite different from 
the small phase lag angle, caused by viscous dissipation, 
between a planet's induced quadrupole moment and the 
perturbing tidal field. A closer black-hole analogy to a 
planet's viscous phase lag may be found by considering 
the tide raised on the hole's horizon by an orbiting moon. 

In this section, we compute the tidal phase shift on the 
horizon for our perturbed Schwarzschild black hole. We 
will discuss in what sense it is and is not analogous to the 
fluid-planet's viscous phase lag. To calculate this phase 
shift, it is convenient to use the Newman-Penrose for- 
malism (see, e.g., chapter 1 of Ref. [13 for a review 
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of the Newman-Penrose formalism) . Appendix [D] sum- 
marizes some details of the Newman-Penrose formalism 
that are relevant for our purpose. 

We consider two approaches to defining the tidal phase 
shift. In Sec. IV Al we define the phase shift only in terms 
of quantities on the horizon (following the method sug- 
gested in Sec. VIIC of Ref. while in Sec. EE wc 
define a phase between the tide raised on the horizon 
and the "retarded" position of the moon (following the 
method used by Hartle in Ref. [l5j'). 



A. Phase of the tidal bulge on the horizon 

For Sec. ^ and Appendix ^| only, we use ingo- 
ing Eddington-Finkelstein coordinates (V", r, 9, <j>) and a 
(+ — — — ) signature of the metric. The Schwarzschild 
metric in these coordinates is 



ds 2 



2M 



dV 2 - 2dVdr 



(46) 

The ingoing Eddington-Finkelstein null time coordinate 
V is related to the Schwarzschild time coordinate t and 
radial coordinate r by the following equation (Eq. (lb) 
of Box 32.2 of Ref. Hfl): 



V = t + r* = t + r + 2Mln \r/2M - 1| 



(47) 



and the Eddington-Finkelstein radial and angular coor- 
dinates {r, 6, <p} are identical to those of Schwarzschild. 

Our slowly orbiting moon deforms the Schwarzschild 
event horizon. By analogy with Newtonian tides, we 
would like to describe the horizon deformation as a per- 
turbation that co-rotates (at a slightly different phase) 
with the tidal field that drives it. But this viewpoint 
inherently envisions the perturbed event horizon as a 
two-dimensional, evolving surface, rather than as a three- 
dimensional, global surface in spacetime. Therefore, be- 
fore we can consider the phases of quantities on the hori- 
zon, we must first specify what we mean by time on the 
horizon. 

Begin by considering the Schwarzschild event hori- 
zon (which is, of course, the three-surface r = 2M). 
There is a preferred way to slice the horizon into a 
single-parameter family of two-surfaces; this preferred 
slicing uses two-surfaces that are orthogonal to the 
Schwarzschild Killing vector d t = d v that is timclike 
at the moon and null on the horizon. Following Har- 
tle , we call this family of two-surfaces the "instanta- 
neous horizon". The instantaneous horizon can be pic- 
tured as an evolving two-surface defined by r = 2M and 
V = const, so that V plays the role of a "time" coordi- 
nate. Throughout this section, we use the terms "hori- 
zon" and "instantaneous horizon" interchangeably unless 
otherwise indicated. 

We now consider how the horizon's perturbation 
evolves with time V. The moon's tidal field, character- 



ized in the LARF by [Eq. 



£■ 



£ 



cxt,20 



£ 



ext,22 
ij 



+ £. 



ext. 2-2 



deforms the otherwise spherical, static horizon. Because 
£.ext,20 j g g .j. a y c anc [ axisymmetric, it cannot contribute to 

the phase shift. For the remaining tidal fields, , 
we shall consider only the (2, 2) mode in detail and the 
result for the (2, —2) mode follows immediately. 

On the horizon, it is the tangential-tangential compo- 
nents of the perturbing tidal field that drive the defor- 
mation (see, e.g., Eq. (6.80) of Ref. ,23); knowledge of 
these components is physicall y eq uivalent to knowledge 
of the Teukolsky function * (see, e.g., Eq. (A7) of 
Ref.jsil). The Teukolsky function is a particular compo- 
nent of the Weyl tensor [Eq. jD5aJ] . 

The horizon deformation is governed by the Newman- 
Penrose equation (Eq. (2.11) of Ref. [T^ l 



(d v - 2e)a 



(i) 



-2(ifi- 



(1) 



(48) 



This is also the "tidal force equation" (equation (6.80) 
of Ref. Here a = — is a Newman-Penrose 

spin coefficient [Eq. l|D3e() ] and Yp-^ is the shear (i.e., the 
rate of change of the shape 5 ) of the instantaneous hori- 
zon. Note that because X and \&o vanish on the unper- 
turbed instantaneous horizon, the spin coefficient e takes 
its Schwarzschild value, which (in our tetrad) is the sur- 
face gravity of the instantaneous horizon gn = (4M) _1 . 

Knowing we can evaluate the horizon shear. Be- 
cause \&q is first order in the perturbation, it may be 
evaluated on the horizon simply by letting r go to 2M. 

Beginning with the (£, to) = (2, 2) metric perturbation 
[Eq. (|12f> . except we now choose the metric signature to 
be (+ — — — )], we compute the perturbed Ricmann 

tensor near the horizon and read off the component ^g 1 ^. 
The result is 



(i) _ • fnftMto V 22 



5 b 3 
= |*^| exp 



y22 e -2ifiV+(8/3)A/Q + Qt]y[ 2^2-j 



(j) - VlV + -MQ - ^ 



-o{M 2 n 2 



(49) 



Here 2X 22 is the spin- weighted spherical harmonic 



^ 2 = W>Mi 



, 2f</ 



With ^ in hand, we can calculate via Eq. I|48|l. 



Recall that the shape of the instantaneous horizon (a two- 
dimensional surface) is completely specified by its intrinsic scalar 
curvature 
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Inserting e = 2gn and u = 2il into Eq. (|48(l yields 



i) 



ico + 2e 



4A/*o 'e 



where 



(l).-2ifi/ 9H 

E«|exp 

>rO(M 2 n 2 ). 

^Horiz 



O(ATfr) 

- fit/ + -m n 



n/g H = 4Mfi. 



7T 

4 



^Horizo 



(50) 



(51) 



T/ie shear X^ 1 ) Zearfs ^q 1 ^ (or equivalently, the perturbing 
tidal field at the horizon) by an angle ^Horizon- [Note 
that the first equality in Eq. ([50"|l appears in Ref. [H| 
Eq. (2.12).] 

The shear is the time derivative of the shape. There 
fore, the shape has a phase 



as 



exp 



2i [<j)-nv 



-Am 



(52) 



In other words, the shear leads the shape by 7r/4. 

The horizon phase shift in Eq. (|51|l follows directly 
from the tidal force equation l|48(l . It is gauge- invariant 
since it only makes reference to gauge-invariant quanti- 
ties measured on the instantaneous horizon. In Ref. j22| 
[Eq. (7.45), Fig. 57, and the surrounding discussion], 
an analogous horizon phase shift <5jjorizon was deduced 
from the tidal force equation for a slowly rotating black 
hole perturbed by a stationary, axisymmetric tidal field — 
physically the same problem as Hartle studied 0] : 



^Horizon = / 9H = 4MQ H = 5 H or 



zion In-.Hff 



(53) 



Here Qh is the horizon angular velocity. 

Although Hartle also used the tidal force equation in 
his calculations, he chose to define the tidal phase shift 
in a different way and made his result gauge-invariant by 
making a connection between the angular positions on 
the horizon and angular positions at infinity through a 
null ray — a choice we will consider in detail in Sec. IV B II 
and apply to our problem in Sec. 

The phase lead (^Horizon i s : m some ways, analogous to 
the phase shift of a tide raised on a non-rotating fluid 
planet. In the latter case, viscous dissipation causes the 
shape of the planet's surface to lag the normal-normal 
component of the perturbing tidal field by a small angle 
<5 v j sc ; somewhat analogously, the horizon shear leads the 
tangential-tangential component of the perturbing tidal 
field. Both phase shifts are small angles associated with 
dissipation (which manifests itself as a secular evolution 
of the energy and angular momentum of the moon's or- 
bit). In the absence of dissipation, there is no phase 
shift. On the other hand, the phase shift (^Horizon is a 
lead angle while S v i sc is a lag angle. Hartle explains this 
difference as a consequence of the teleological nature of 



the horizon ^f|. Also as Hartle observed, when the an- 
gular velocity fl is not small compared with 1/M, the 
deformation of the horizon cannot be described in terms 
of a phase shift [Til ]. 



B. Phase shift between the tidal bulge and the 
moon 

As an alternative to the above way of defining the tidal 
phase shift, one can define it as the angle between the 
tidal bulge on the horizon and the location of the moon 
in its orbit. Hartle used this approach when he computed 
the tidal lead on a rotating hole perturbed by a stationary 
moon First, we will briefly summarize the aspects 
of Hartle's analysis which are relevant to our purpose. 
Then, we will apply his method to a slowly rotating moon 
around an otherwise Schwarzschild black hole. 



1. Tidal phase shift between a rotating horizon and 
stationary moon 

In Ref. ^5|, Hartle considers the problem of a distant, 
stationary moon perturbing a slowly rotating black hole. 
The Kerr metric can be written as 



ds 2 = 



2Mr\ 
1 - dV 2 



AaM r sin 2 9 - , - 
2dVdr H dVdcj) 



-2a sin 2 9drd<p - Ed9 2 

, 2a 2 Mr sin 2 9 



(54) 



Here £ = r 2 + a 2 cos 2 9. The coordinates V and (f> are 
related to the usual Boyer-Lindquist coordinates t and <p 
by 



dt 



dV - 



-dr 



A 



-dr 



(55) 



where A = r 2 — 2Mr + a 2 . When a = 0, Eq. l(5*4*|) reduces 
to the Schwarzschild metric in Eddington-Finkelstein co- 
ordinates [Eq. (|4^|) ]. 

The event horizon is the surface r = r + = M + 



VM 2 ^ 



Just as in the Schwarzschild case consid- 



ered above, the event horizon can be sliced into a single- 
parameter family of two-dimensional surfaces using the 
Killing vector dy which is timelike at infinity and null on 
the horizon. This family of surfaces is the instantaneous 
horizon. 

The distant moon raises a tidal bulge on the central 
hole's instantaneous horizon. In the limit that the moon 
is far away, the change in the horizon's shape (or equiv- 
alently, the change in the scalar curvature of the in- 
stantaneous horizon), is purely quadrupolar. 
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The deformation is driven by the transverse-transverse 
component of the tidal field at the horizon, which is phys- 
ically equivalent to the Tcukolsky function, a particular 
component of the Riemann tensor ^>o [Eq. l|D5a(l ]. This 
component vanishes in the unperturbed Kerr spacetime 
[Eq. (|D12a|) ]. and the first order correction vffg 1 ^ has the 
form 



(i) 



S' m (r) 2 Y 



Ira i 



(56) 



the instantaneous horizon at angular position 6 
= a/b-a/2M. 



moon 



(60) 



The tidal bulge therefore leads the moon's position by an 
amount 



null map 



6 H 

moon 



bulge 



M 2 n H 



(61) 



where 2 Y em is a spin-weight-2 spherical harmonic. Be- 
cause the perturbation is purely quadrupolar, we need 
only consider the case I = 2, m = 2 here, although Hartlc 
considers the generic case. Hartle uses Teukolsky's solu- 
tion [35| for the stationary radial functions S tm due to the 
£-polc perturbation caused by a distant, stationary point 
particle with mass fi. Furthermore, while Hartle treats 
the case of a moon at any location (9, </>), for concreteness 
we specify the moon's position as (0,<fr) = (w/2,0). On 
the horizon, the Teukolsky function turns out to have the 
value (combining Eqs. (4.30), (4.31), (4.15), and (4.18) 
of Ref. 0) 



2i 



(!) H iflMflH . 4 , 

= vi^ sin U ,cxp 



2MQ.H 



(57) 



The tidal field deforms the instantaneous horizon, 
changing it's shape and thus its two-dimensional scalar 
curvature S%. Hartle computes the quadrupolar correc- 
tion to the scalar curvature, ^ , ^ 1 - ) ' f=2 ' H , of the instanta- 
neous horizon [Eq. ljD9|) ]. His result is (Eqs. (4.26)-(4.27) 
of Ref. 0) 



(l),^=2,H 



OC COS 
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H 



+0 i^pj + O (M^ H ) . (58) 



Here we have used the relation (valid for small a/M) that 
a = 4M 2 flH, with being the angular velocity of the 
hole. For simplicity, one can then take the limit b — > 00. 

Before continuing, we should remark that Hartle's pre- 
scription for constructing S^ ull m can be described with- 
out reference to the moon's position. Begin by comput- 
ing the angular location of the tidal bulge on the hori- 
zon. Next, ingoing, zero-angular-momentum null rays 
from infinity define lines of constant angle, so that there 
is a one-to-one correspondence between angular positions 
on the horizon and angular positions at infinity. The 
angular position at infinity of the tidal bulge can thus 
be computed. Finally, perform the calculation again, 
but this time perturb a non-rotating spacetime; in this 
case, there will be no tidal friction. Because the Kerr 
and Schwarzschild spacetimes are asymptotically identi- 
cal, one can unambiguously compare the angular position 
of the tidal bulge in the presence and in the absence of 
tidal friction: Jmhii = ^,,1^ ^,,1^ nn fr^+;^ n - This 



map 



pro- 



of A H 

Ui "null map 



^bulgc ^bulge, no friction 

is equivalent to the previous definition of S^ uii 
vided that b — > 00. 

However, this alternative formulation 
breaks down when the moon, not the horizon, rotates 
The rotation is then described by Q, which is a parame- 
ter of the perturbation, not of the background spacetime. 
To eliminate tidal friction, one must let Q — > 0, which 
eliminates the perturbation 7 . Because of this failure, we 
prefer to consider Hartle's phase shift as a comparison of 
the position of the tidal bulge with the position of the 



Instead of measuring the angle between the shear a 
and the tidal field "to on the horizon, Hartle defines 
his phase lead as the angle between the shape and 
the moon's angular position. To make this definition 
gauge- invariant, Hartle chooses ingoing, zero-angular- 
monentum, null geodesies to be "lines of constant angle." 
He then compares the angular position of the horizon 
tidal bulge, 

to the angular position of the moon on the horizon. 

Consider stationary moon in the equatorial plane at 
(large) radius r = b and at angular position — 0. An 
ingoing null ray, originating from the moon, intersects 



2. Tidal phase shift between a non-rotating horizon and 
rotating moon 

A similar analysis can be applied to our system, in 
which a distant moon in a slow, circular orbit raises a 
tide on a non-rotating black hole. 



° Note that there is a sign error in Hartle's analysis. Hartle in- 
correctly states that the ingoing null ray intersects the horizon 
at +a/2M + 0(a/b), not -a/2M + 0{a/b). Had wc also made 
this error, there would be a coefficient of 20/3 instead of 8/3 in 
Eq. l6ll . 

7 Even if f! — > resulted in a non-zero perturbation, it is unclear 
how to distinguish such a perturbation from a small change in 
the coordinates of the background spacetime. 
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The moon orbits the central black hole along the world 
line specified by Eq. ©. In other words, the moon has a 
phase given by 

</>moon(V0 = 0* = n(V - 6*). (62) 

This must be compared with the location of the bulge 
on the hole's future horizon. Equation (|52fl for [or, 
alternatively, inserting Eq. (|49|l into Eq. l|D9(l ] shows that 
the tip of the tidal bulge has a phase given by 

g 

<?!>buigo = fiv + -Am. (63) 

As time V passes, this bulge rotates around and around 
the horizon, with the same angular velocity il as the 
moon that raises the tide. 

Following Hartlc, we compare the angular location of 
the tidal bulge, <^buigc(V), with the angular location of 
the moon, </> m0 on(^), using ingoing, zero angular momen- 
tum (ZAM) null rays to provide the connection between 
<f> at the moon's orbit and <f> on the horizon. In the ingo- 
ing Eddington-Finkelstein coordinates that we are using, 
these ZAM rays have a very simple form: 

{V", 6, <fr} = constant , r decreases from b to 2M . (64) 

Since V, 8 and <j> are all constant along these rays, they 
give us a one-to-one map of events {V, r = b,6, </>} at the 
moon's orbital radius to events {V, r = 2M, 8, 0} on the 
horizon that have identically the same V, 8, and <f>. With 
the aid of this map, we conclude that the angle by which 
the horizon bulge lags the moon's position is 

8 

<5null map = 0bulge(V) - (/> m oon{V) = -Mil + fib*. (65) 

[Eqs. (|62|l and Again, the phase shift is a phase 

lead, not a phase lag, due to the tclcological nature of 
the horizon. 

In addition to the teleological phase shift of order Mil, 
5 nu \i map contains a much larger term of magnitude fib* ; 
this term reflects the choice to use an ingoing-null-ray 
mapping between the moon and the horizon. A similar 
term appears in Hartle's calculation [Eq. I|61(l ]. but in 
Hartlc's system the term is much smaller than the tele- 
ological phase shift size (specifically, smaller by a factor 
of M/b), whereas fib* > MCI. 

One could avoid this problem by defining the phase 
shift to include only terms of order Mil and Mflji. With 
this definition, the remaining tidal phase leads are the 
same: (8/3)Mf2, as one would expect, given that a there 
should be no tidal shift at all if the moon were to rotate 
at the hole's angular velocity, i.e., if fl = fljj. 

We prefer, however, to define the tidal lead angle in as 
the angle ^Horizon = 4Mf2 by which the horizon shear 
leads the horizon tidal field. This angle, in contrast 
to <5 nu ii map , is defined in terms of an "instantaneous" 
(spacelike) connection between the moon and the hori- 
zon, i.e., by the near zone mapping of the moon's posi- 
tion to the horizon tidal field's [\&q 's] maximum. [Had 



the moon been in the radiation zone (b ^> \/2ir), one 
would have expected the connection to be lightlike.] 



VI. CONCLUSION 

For our simple system of a Schwarzschild black hole 
and circularly orbiting moon, we have found that the 
time-dependent part of the moon's tidal field induces a 
quadrupole moment that is unambiguous. The static 
induced quadrupole moment was found to be zero in 
the Regge- Wheeler gauge, but it is ambiguous in gen- 
eral. The ambiguity of the static induced quadrupole 
moment leads to an ambiguity in the phase of the in- 
duced quadrupole moment in the LARF; however, the 
tidal bulge on the horizon still has a well defined phase 
shift with respect to the orbiting moon. Because of the 
ambiguity of the induced quadrupole moment and the 
LARF phase shift, we conclude that the polarizability 
and phase shift are not suitable for constructing a body- 
independent description of tidal coupling in EMRIs. 

However, this conclusion does not eliminate the possi- 
bility of developing a body-independent language to de- 
scribe tidal coupling, including cases where the central 
body is a black hole. It might be possible, for instance, 
to define a new set of induced "dissipative multipole mo- 
ments" for the central body — i.e. moments that vanish 
in the absence of tidal friction. Such dissipative moments 
would still be linear in the perturbing tidal field, so one 
could still define a polarizability. Also, by ignoring any 
non-dissipative tidal coupling, the phase shift might no 
longer contain additional information. Even if such an 
extension does not prove feasible, tidal coupling can still 
be described in the more conventional (but still body- 
independent) language of energy and angular momentum 
transfer between the moon and the central body. 

Other future work could include generalizing our anal- 
ysis to spinning black holes, treating noncircular, non- 
equatorial orbits, and (most importantly) studying how 
information about tidal coupling in EMRIs can be ex- 
tracted from the gravitational waves detected by LISA. 
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APPENDIX A: SYMMETRIC TRACE-FREE 
TENSOR NOTATION FOR SPHERICAL 
HARMONICS 

The scalar spherical harmonics Y im (9, <f>) can be writ- 
ten in terms of of rank-^ symmetric trace-free (STF) ten- 
sors [23|. The spherical harmonics Y 2m (9,(j)) that have 
been used in this paper are 




±2i4 



3cos26>) . 



(Ala) 



(Alb) 



They can be written in terms of rank-2 STF tensors as 
(Eq. (2.11) of Ref. [13) 

Y^(0 t <f>) = 9@ n n i n j , (A2) 

where n l = x % jr and £^ 2m are the STF tensors given by 
(Eq. (2.12) of Ref. gj)' 
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2±2 




(A3) 



(A4) 



APPENDIX B: TIME-DEPENDENT 
PERTURBATION EQUATIONS 

In Regge- Wheeler gauge, the metric perturbation for 
a given even-parity (£, to, lu) mode depends on the three 
radial functions H, H\, and K. In this appendix, we 
introduce the Zerilli function Z and the Regge- Wheeler 
function X and describe how we obtain the radial func- 
tions from them. The description here will hold for a 
general (£,m,u>), while the results derived in Sec. I|III|) 
rely on the special case when (£, to, lo) = (2, 2, 2S7). 

The original Zcrilli's master function is defined implic- 
itly through its relation with the two radial functions Hi 



and K [Eqs. (13) and (14) of Ref. 
placed by Z]: 



with r l 



M 



H\ = —ILO 



Xr 2 - 3XMr - 3M 2 



K 



(r - 2M)(Xr + 3M) 
A(A + l)r 2 + 3XMr + 6M 2 
r 2 (Ar + 3M) 



dZ 

Z — iujr——, 



dr 
dZ 
dr* 



(e) ^ 

(Bl) 
(B2) 



where 



X= -(£-!)(£ + 2). 



Using the algebraic relationship (Eq. (10) of Ref. |2(j) 



H = 



i(X + l)M 



Hi 



x M 

X + — 

r 



M 



'/r 2 + bj 2 r 2 
1 - 2M/r 



K 



one can obtain H in terms of the Zerilli function 



H = 



in which 



si 



2M-r r 2 (3M + Xr) 2 



Z + s 2 — , (B3) 
dr 



si = 9A/ 2 (7\/ + Ar) + A 2 r 2 [3M+(A + l)r], 
-MI 2 - 3XMr + Xr 2 



S2 



r(3A/ + Xr) 



The Zerilli function obeys the wave equation (Eqs. (18) 
and (19) of Ref. 0): 



dr* 2 



+ uj 2 - V(r) 



Z = 0, 



(B4) 



in which the potential term is given by 
2(r-2M) 



V(r) 



r 4 (Xr + 3M) 2 



A 2 (A + l)r 3 + 
3A 2 Mr 2 + 9AM 2 r + 9M 3 



The odd-parity master function, the Regge- Wheeler 
function, is defined in Eq. (23) of Ref. [24| (and is called 
Q in Regge and Wheeler's notation). It obeys the differ- 
ential equation (Eq. (7) and of Ref. [2(|): 



dr* 2 



_L) 6M \ 



r 3 J 



X = 0. 



(B5) 

The connection between the Regge- Wheeler and Zerilli 
functions was first found by Chandrasekhar and is listed, 
e.g., in Eq. (152) of Ch. 4 of Ref. 0: 



A(A+ 1) -3iMw 



S 



9M 2 (r-2M) 
A(A + 1)+ r 2 (Ar + 3M) 
2M\ dX 
dr 



3M 1 



(B6) 



This completes our metric reconstruction scheme from 
the Regge- Wheeler function. We are now ready to evalu- 
ate the radial metric perturbation functions H, Hi, and 
K for the (£,m) = (2, 2) mode of the perturbations. Ex- 
panding X H [given in Eq. (|16|l in powers of f = r/2M to 
first order in e = 2Mlo, we obtain 



X H =A 
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-3 



l + T2 l£ F + ^ 

n— 5 



IS 



+ 0(e 2 ) 



(B7) 



Here A is an overall scaling factor (Eq. 1)20(1 in Scc. lIII A| . 
While the summation can be rewritten as a closed-form 
expression, we prefer to stay in the series notation, since 
our interest is in reading various powers of r in the re- 
sulting first-order tidal field. Equation l|B7|) is the value 
of the Regge- Wheeler function in the LARF; inserting it 
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into Eq. (|B6(1 yields the expression for Z in the LARF 
[We shall suppress "0(e 2 )" hereafter]: 



1 + 



A 



3r 2 



63 



9r 21 

4~ ~~ 16 ~~ 64 + 256f 
945 - 236ie 8505 + 15436ie 



5120f 2 61440f 3 

Inserting Z into Eqs. (jEJ, (JB2J) and (jHSJl yields A' 
and H . Expanded in powers of f and to first order in e, 
these radial functions are given by 



H = 
H x = 

K = 



A_ - 

M . 
iAe 
4M 

A_ r 

M 



(3 + 4ie)(f 2 -f 
- 8f 3 - 2r 2 + 



4f 



f- 1 



r 2 



(3 + 4ie)(r 2 



l£ 

8r* 



1 

2 

70f5 



1 


10r 3 




+ 1 




- 3 + 


r 


ie 


+ 10f 3 




+ 







3/c 



20f 4 



Of 



+ o(f- b ) 



+ 0(f~') (B8) 



APPENDIX C: TIME-INDEPENDENT 
PERTURBATION EQUATIONS 

As is evident from the time-dependent perturbation 
theory, as co — ► 0, i?i goes to zero. In the static case, 
then, there are only two radial functions, H and K 
(where the hats signify that they represent the time- 
independent perturbations). Specializing to the axisym- 
metric case, the metric perturbation is 



l ab 



Diag (l-2M/r)H, 



H 



r 2 K, r 2 K sin 2 I 



1 - 2M/r 
x Y m (8. 



The linearized Einstein equations governing H and K are 
given in Eqs. (9d) and (9e) of Ref. [U with H i = and 
u> = (k = in Edclstcin and Vishvcshwara's notation): 



dK 
dr 

2M dK 



dH 

dr 



1 



2M 

2 



r 
2M 
r 



2M^ 1 

d 2 H 2 dH 
r dr 



1 



dr 2 



H, 



1) 



(CI) 

H. 

(C2) 



Eliminating dK / dr from these two equations, we can 
then write a single second-order differential equation for 
H in terms of the variable z = r/M—1 (same as Eq. (4.9) 
of Ref. [3): 



(1-z 2 ) 



d 2 H 
~dz^ 



2z- 



dH 

dz 



1(1 + 1) 



1-z 2 



H = 0. 



This takes a form of the associated Legendre differential 
equation. The general solution for H is therefore 

a t P 2 (r/M - 1) + l3 t Q 2 (r/M - 1). (C3) 



With the general solution for H, we can integrate 
Eq. iJHII or (dl to find K. For I = 2, we have 

K 2 (r) = [a 2 P*(r/M-l)+f3 2 Ql(r/M-l)] 

x ^=H7! + ffi(r) ' (C4) 



APPENDIX D: NEWMAN-PENROSE 
FORMALISM 

In this appendix, we summarize some equations of 
the Newman-Penrose formalism for our choice of tetrad. 
In this Appendix and in Sec. ^ only, we use ingo- 
ing Eddington-Finkelstein coordinates (V, r, 8, <j>) and a 
(+ — — — ) signature of the metric. 



1. Newman-Penrose quantities for Schwarzschild 

We adopt the Hartle-Hawking null tetrad, which is 
given by Eqs. (4.2) of Ref. together with the nor- 
malization conditions l^n^ = 1 and rn^fn^ = — 1. The 
tetrad vectors have components [using the notation e" = 

(eV^e*)] 



n" = (0,-1,0,0) 
m" = fo, 0, — - 

'0,0, 



V2r ' y/2r sin (9 
1 i 



v^r' v^rsinfl 



(Dla) 
(Dlb) 
(Die) 

(Did) 



Note that throughout this Appendix, an overbar denotes 
complex conjugation. 

From these basis vectors, we define the direction 
derivatives 

D = A = n^, 6 = m^d^ and 5 = m"d M . (D2) 

Our conventions for the spin coefficients follow Ref. C3 
[specifically, Eqs. (2.2) and (2.3)]. The spin coefficients 
arc defined by 

(D3a) 
(D3b) 
(D3c) 
(D3d) 
(D3e) 
(D3f) 

(D3g) 



H 



K = 






n = 






P = 


t fl . v m fl m' / 




fj, = 


—n^. v rh ,1 m v 




a = 






A = 


—n^ v fn} l fh v 




e = 




- m KV m*t) . 


a = 


X - {l K „n»rh v 


— m fl - v fh fJ, ifi l/ ) 


P = 


\ (i^m v 


— m ll - v fnf i m v ') 
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The spin coefficients for the Schwarzschild spacctime 
are 



a = X = v = t 

M 

2T 2 
r-2M 



7T = 7 = 



p 
II 



2r 2 



1 

r 

-/? = 



2\/2rtan0' 



(D4a) 
(D4b) 

(D4c) 

(D4d) 

(D4e) 



Because we are only interested in vacuum regions of 
spacctime, the Riemann and Weyl tensors are inter- 
changeable. The Weyl components are defined in vacuum 
by 



*0 = 




(D5a) 


*1 = 




(D5b) 


*2 = 


—^R a /3yS (^n^Pn 6 + l a rv a m~ 1 m i ) 


(D5c) 


*3 = 


-R a 1 d a n l3 m'n s 


(D5d) 


* 4 = 


i > a — 3 -y — 6 

—RaB-ysn m H n'm . 


(D5e) 



Their values for the Schwarzschild spacctime are 

*o = *i = *3 = *4 = 0, 
M 

The Ricci identities are 



(D6a) 
(D6b) 

(D7) 



Inserting the null tetrad vectors for e a and projecting 
along the tetrad yields the Ricci identities in Ncwman- 
Penrose notation. One of these equations is, in our tetrad 
and evaluated on the horizon. 



Da™ - 2eaW 



2ea (1) = * 



(i) 



(D8) 



(Note that we have used the fact that a and \&o vanish for 
Schwarzschild.) This is the tidal force equation; it relates 

^Pq 1 ', which is physically equivalent to the tangential- 
tangentail component of the perturbing tidal field on the 
horizon, to a, which is physically equivalent to the shear 
of the instantaneous horizon. 

The shape of the perturbed instantaneous horizon 
is determined by its two-dimensional scalar curvature 
3t + &W where 3$ is the curvature of the unperturbed 



horizon. According to the tidal force equation (|D8(I . ^ 
drives the shear, which is the "rate of change of the 
shape" of the horizon as measured by fiducial observers 
on the horizon J^. Thus, it is not surprising that 

can be computed directly from Vl/g 1 . Hartle [l^ has de- 
rived [T5) the explicit formula, a consequence of Gauss' 



relation |36| , in the Ncwman-Penrose formalism with the 
present choice of coordinates and tetrad: 



-41m 



(S + 2n- 2a)(S + tt - 4a) + 2eA 



Lo(iu> + 2e) 



(i) 





(D9) 

where w is the frequency of the perturbation. When a 
Schwarzschild black hole is perturbed by a distant moon 
in a slow, circular, orbit with angular velocity f2, then 
lj = 2fl. 



2. Newman-Penrose quantities for Kerr 



Finally, to facilitate our comparison to Hartle's results, 
we here list the relevant Newman-Penrose quantities for 
the Kerr spacetime [Eq. (|54() ] using Hartle's choice 0| 
of coordinates and tetrad. In the limit a = 0, Hartle's 
tetrad and spin coefficients reduce to those listed in the 
previous subsection. 

The null tetrad vectors [using the notation e p = 

(e 9 ,e r ,e e ,e^)} are 



2Mr + a 2 a 

" 1 ; 



2 (r 2 + a 2 ) ' ' r 2 + a 2 
2{a 2 + r 2 ) 



(DlOa) 



2r 2 + a 2 + a 2 cos 29 



0,0 



1 2 cos 29 



2 {a 2 + 2r 2 + a 2 cos 29) 
—asm.9 



-HI 



\[2{ir + a cos 9) 

—asin.9 

■—= m M 

V2 (—ir + a cos 9) 

asinfl (r 2 - 2Mr + a 2 ) 



(DlOb) 



mr 



0. 



2\f2 (r 2 + a 2 ) (-ir + a cos 9) ' 
1 (ir + a cos 9) esc 9 



V2 (r + ia cos ( 



V2{r 2 +a 2 



(DIOc) 



0,- 



asinfl (r 2 -2Mr + a 2 ) 
2V2 (r 2 + a 2 ) (ir + a cos 9) ' 
1 (— ir + a cos 9) csc( 



V2( 



r — ia cos t 



V2(r 2 + a 2 ) 



(DIOd) 

Then, one can compute the spin coefficients for this 
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tetrad from Eqs. i|D3a|) (|D3i 



The directional derivatives arc then given byEq. (lD2l . 



o 



K 

A 



(7 = 

0(a 2 ) 
0(a 2 ) 

-i(2M + r) sin 8a 

i(4M + r)sin 9a 

—i cos 8a „, 9% 
+ 0(a 2 ) 



2r 2 
M r 2 -a 2 



M 



P = 



r-2M i(r -2M) cos 9a 
2~r~ 2 2~?> 

-0(a 2 ) 

- + 0(a 2 ) 



(Dlla) 
(Dllb) 
(Dllc) 

(Dlld) 
(Dlle) 
(Dllf) 



(Dllh) 
(Dili) 



cot 9 i[-3M + (2M + 2r) cos 29]a 



4V2r 3 sin6> 



2V2r 

+0(a 2 ) (Dllj) 
cot0 i[M + r + (r-M) cos 26]a 



2\f2r 
+0(a 2 ). 



4V2r 3 sin 9 



(Dllk) 



Using the Kerr metric [Eq. 1|54|) ] and Hartle's choice 
for the tetrad [Eqs. l|D10ajl -l |D10d|) ]. one can compute 
the curvature for Kerr and read off the curvature scalars 
via Eqs. (|D5a|) - (|D5e|) : 



^4 



*i = 



M 



(r — ia cos 9) 
3iaAI s'm9 

\2{r — ia cos ( 
3ia 2 M sin 2 9 
(ir + acos9) 5 ' 



(D12a) 
(D12b) 

(D12c) 
(D12d) 



The tidal force equation (|D8j) relates to er (1) . The 
correction to the scalar curvature of the horizon, , is 
given by Eq. (ID9I) . 



For a stationary moon perturbing a slowly rotating 
Kerr black hole, the frequency of the perturbation is lu = 
-2Q.H = -8M 2 n H . 
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